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The  Complete  Electromagnetic  Fields  in  the 
Focal  Region  of  a Paraboloidal  Reflector 


1.  INTROWCTION 

In  designing  optimum  feed  systems  for  Cassegrain  reflector  systems,  it  is 
highly  desirable  to  have  an  accurate  picture  of  the  electromagnetic  fields  in  the 
reflector  focal  region.  In  order  to  study  these  fields  we  have  considered  the  case 
of  a plane  wave  incident  upon  a large  reflector,  as  shown  in  Figure  1,  and  have 
used  the  physical  optics  approximation  to  calculate  the  complete  electromagnetic 
field  distribution  produced  in  the  vicinity  of  the  reflector  focus. 

2.  ' THEORETICAL  BACKGROOD 

Let  us  consider  a plane  wave  with  electric  and  magnetic  fields* 


Ej  ■ y exp  [ i (wt  + kz) ) 
X exp  [ i(wt  + kz)  ) 


(Received  for  publication  29  September  1976) 

*In  Eq.  (1)  X,  y,  and  z are  unit  vectors  along  x,  y,  and  z.  Also  k is  the  wave- 
number  * 2»/A,  where  ^ is  the  signal  wavelength. 


incident  from  the  right  upon  the  reflector  in  Figure  1.  If  we  assume  that  the 
reflector  surface  is  described  by  the  arbitrary  function 


z = f(x,y) 


(2 


it  can  be  shown^  that  the  magnetic  field  scattered  by  the  reflector  is  given,  in  the 
physical  optics  approximation,  by 


HS 
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dy 


* ikz 
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where  R is  the  distance  from  a source  point  (x,  y,  z)  on  the  reflector  to  the  field 
point  (Xp,  y^,  z^),  is  the  projection  of  the  reflector  surface  onto  the  x - y plane 
and  x.y,  and  z are  unit  vectors. 
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Figure  1.  Reflector  Geometry 


1.  Silver,  S.  (1965)  Microwave  Antenna  Theory  and  Design.  Dover  (New  York). 


If  we  assume  that  the  reflector  surface  is  a parabola  with  a focus  at  z > F, 
then  Eq.  (2)  becomes 


4F 

and  the  projection  of  the  reflector  onto  the  x-y  plane  is  a circle  satisfying  the 
equation 


(4) 


2 2 

•X.  + y 
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where  D is  the  diameter  of  the  reflector.  If  we  now  use  Eqs.  (4)  and  (5)  in  (3)  we 
obtain,  after  some  manipulation 
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where 


«(x,y)  ^)exp[-ik  j. 

R = i<x-x^)2  + (y-y^)2  + . 


The  electric  field  distribution  can  be  obtained  by  employing  the  Maxwell  equation 


VX  - iuJC^Eg 


The  result  for  Eg  is 


(7) 
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where 


2 2 y y 


Equations  (6)  and  (8)  are  the  formal  expressions  for  the  complete  electro- 


magnetic  fields  in  the  physical  optics  approximation.  They  represent  a quite  good 
approximation  for  the  entire  region  z > (x^  + y^).  which  is  of  interest  to  us. 

Of  course,  they  are  inaccurate  for  z < 0;  over  part  of  that  region  the  geometrical 
theory  of  diffraction  must  be  employed. 

By  observing  Eqs.  (6)  and  (8)  it  is  clear  that  the  scattered  electric  and 
magnetic  fields  possess  certain  symmetry  properties.  These  are  (for  a fixed  z^) 
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Because  of  the  aforementioned  symmetry  properties  we  have  calculated 
E and  H only  for  positive  values  of  x^  and  y^;  the  values  for  negative  x^,  y^  follow 
immediately  from  Eqs.  (9)  and  (10). 


3.  RESULTS 

We  have  developed  a computer  program  to  calculate  the  field  components 
given  by  Eqs.  (6)  and  (8).  As  an  example  of  typical  results  of  our  program,  we 
have  studied  a reflector  such  that 


X . A. 

D 3 


T ' 60 


where  A is  the  signal  wavelength,  and  have  calculated  the  field  distribution  in  the 
planes  = 0.  95F,  0.  967F,  0.  983 F,  and  1.  OF.  In  Figures  2 and  3 we  show*  the 
amplitude  and  phase  of  the  electric  and  magnetic  fields  in  the  plane  z^  > 0.  95F. 
The  fields  shown  are  those  along  the  line  y^  ‘ 0,  for  differing  values  of  Xg.  In 
Figures  4 and  5 we  show  the  fields  along  the  line  Xg  • 0 for  differing  values  of  yg. 


*In  all  the  results  of  Figures  2 to  11  we  have  assumed  Hg  > 2z. 


( 


Figure  2.  Magnitude  of  the  Fields  Along 
the  yo  = 0 Axis  for  z = 0.  95F 


Figure  4.  Magnitude  of  the  Fields  Along 
the  Xq  = 0 Axis  for  Zq  = 0.  95F 


Figure  3.  Phase  of  the  Fields  Along  the 
yg  = 0 Axis  for  z^  = 0.  95F 


Figure  5.  Phase  of  the  Fields  Along  the 
Xq  = 0 Axis  for  Zq  = 0.  95F 
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We  note,  upon  comparing  Figures  2 and  4.  that  if  the  results  are  known  along  the 


y =0  aucis  we  can  immediately  obtain  those  along  the  axis  ■ 0 by  replacing 
iH^l  by  lE  1 /Z^.  IhJ  by  IeJ  /Z^.  and  Ie^I  /Z^  by  Ih^I  . where  Z^  is  the 
impedance  of  vacuum.  Because  of  this  duality,  in  the  remaining  figures  we  wUl 
only  show  results  along  either  the  x^  * 0 or  the  y^  ■ 0 axis.  In  Figures  6 and  7 
we  show  the  fields  in  the  plane  ■ 0.  967F;  in  Figures  8 and  9 we  show  the  results 
in  the  plane  z^  * 0.  983F.  and  finally,  in  Figures  10  and  11  we  show  the  results 
in  the  focal  plane.  The  components  E^^  and  are  now  shown  bee.  use  they 
are  both  zero. 

There  are  several  obseirvations  which  should  be  made  regarding  our  results: 

(1)  The  results  of  Figure  10  for  the  transverse  fields  in  the  focal  plane 
agree  with  those  calculated  earlier  by  Minnett  et  al.  ^ (see  Figure  11  of  Minett's 
paper,  for  our  case  74°). 

(2)  The  cross -polarized  fields  Ie^I/Z^  and  1h^|  are  generaUy  of  the  same 
order  as  |H^I  and  |Ey|/Z^.  except  very  near  to  the  z^  axis.  This  is  true  even 
in  the  focal  plane,  and  even  holds  within  the  focal  spot  (that  is.  we  call  the  trans- 
verse dimension  of  the  first  null  in  Figure  10  the  focal  spot  size,  and  this  is  of 
order  Xf/D)  as  can  be  seen  from  Figure  10.  where  1e^|  /Z^  is  small  near  the 
center  of  the  focal  spot  (yo  near  zero)  but  is  large  near  the  outer  edge  (y^~  0.  SX). 

(3)  In  Figures  2 to  11  we  have  shown  the  fields  on  the  x and  y axes  where 

= Hy  = 0.  This  should  not  imply  that  Ej^  = Hy  = 0 off  these  axes.  In  Figure  12 
we  show  the  field  distribution  along  an  axis  (see  Figure  13)  oriented  45°  relative 
to  the  X axis  in  the  plane  Zq.  Note  that  both  and  Hy  are  nonzero,  although  they 
are  considerably  smaller  than  the  other  field  components. 


2.  Minnett.  H.  C.  and  Thomas.  B.  (1968)  Fields  in  the  image  plane  of  symmetrical 
focusing  reflectors.  Proc.  IEEE,  115:14 19- 1430. 
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Appendix  A 

In  this  appendix  we  present  a Fortran  listing  of  the  computer  program  used 
to  calculate  E and  H.  Note  that  the  quantities  printed  out  for  E are  actually  the 
electric  field  normalized  by  rather  than  E.  The  inputs  to  the  program  are: 

D = diameter  of  reflector, 

F = focal  length, 

K = wavenumber  = 2x/X  , 

XO,YO,  ZO  = coordinates  of  observation  point, 

XTOL  = YTOL  usually  set  to  lO"^. 


PR06RAH  NAINI INPUT t OUTPUT) 
IHPLICIT  CONPLEX(Q) 
COHHON/ONE/XO,  ONEOEFt  0NE04F, 
COMMON/QXI/QXINTEG (6) 
COPHON/XTOL/XTOL 
COPNON/VTOL/YTOL 
REAL  K 

NAHELIST/XINPUT/  0«  Ft  K*  Z0«  XOt 
READ  XINPUT 
OCONST  s (1.0*  0.0)  / 

0.5  / F 
0.25  / F 
(0  / 2.0)**2 
2.0 


20.  K,  TO.  OOV2SQ*  OCONST 


TO*  XTOL*  YTOL 


CHPLXIO.O*  K) 

0NF02F 
ONEOAF 
OOV2SO 
XB  » 0 / 

XA  s -XB 

CALL  XINTEGfXB*  XA) 

PRINT  FOO*  QXINTEG 
STOP 

600  F09NAT(6H  HX  s .1P2E12.5/6H  HY  * «2E12.5/6H  HZ  > .2E12.5/ 
1 6H  EX  « (2E12.5/6H  EV  = .2E12.5/6H  E2  > *2E12.5) 

END 


SUBROUTINE  XINTECIXBt  XA) 

INPLICIT  CONPLEXCQI 
CC**PON/OXI/OXIKrEG  16) 

C0PN0N/QVI/0YINTE6(6I 

COHNON/XTOL/TOL 

DIMENSION  QTH0I6)«  0F0UR(6>.  OENOSIbIt  QT0TAL(6) 

H s (XB  - XA)  / 2.0 
N a 1 

CALL  YINTE6IH»XA) 

DO  10 

QTNC(J)  a (0.0.  0.0) 

10  OEOUR(J)  a QyINTEG(J) 

CALL  YINTEG(XA) 

00  20  Jal,6 

20  QENOS(J)  a OYINTEG(J) 

CALL  YINTEG(XB) 

DO  30  Jal.6 

QENDS(J)  a QENDS(J)  ♦ OYIMTEG(J) 

30  QTCTAL(J)  a (QENOSIJ)  * 4.0»QFQUR(J)  ) • H / 3.0 
40  DO  50  Jal.6 
50  QXINTEG(J)  a QTOTAL(J) 

Y a H a H / 2.0 
N a 2 * N 

DO  60  Jal.6 

QTMO(J)  a QTNOfJ)  * OFOUR(J) 

60  QFOUR(J)  a (0.0.  0.0) 

1 = 0 

70  I a I ♦ 1 

CALL  YINTEG(Y^XA) 

DO  80  Jal.6 

80  QFOUR(J)  a QFOUR(J)  *■  OYILTEG(J) 

Y a Y ♦ H « H 

IF(I  .LT.  K)  GO  TO  70 
IFLAG  a 0 
PRINT  * 

DO  90  Jal.6 

QTOTAL(J)  a (QENOSIJ)  t B.C^OTMOIJ)  * 4.0*QFOURIJ)  ) 
QO<:NOH  a OTOTALIJ) 

IF (CABSIQOENOH)  .LT.  TOL)  QOENOF  = CNPLXITOL.  0.0) 
IFICABSI  lOXINTEGIJ)  - OTOTALIJ)  )/OOENON)  .GT.  TOL) 
PRINT  «.  CABS  I OTOTALIJ)  ) 

90  CONTINUE 

IFIIFLAG  .EQ.  1)  GO  TO  40 
DO  100  Jal.6 
100  QXINTEGIJ)  a OTOTALIJ) 

RETURN 

END 


♦ H / 3.0 
IFLAG  a 1 


SUBRO^ITINE  VINTEGIXI 
IMPLICIT  COMPLEX  IQ) 

COPMOM/ONE/XO,  ONEOEF,  0NE04F,  Z0«  K,  YOt  OOV2SQ*  OCONST 
C0MM0N/0Q0/0INTI6) 

C0MM0N/QVI/QYINTE6I6) 

COMMON/XXX/XOIFFt  X0IFF2«  XSOt  EX 

COMMON/YTOL/TOL 

REAL  K 

DIMENSION  QTN0I6)(  QF0URI6lt  0EN0SI6) t QT0TALI6) 

EX  « X 

XSO  * X**2 

XOIFF  * XO  - X 

XOIFF2  » XOIFF**2 

YB  s 6AMX  X SORTIOOV2SO  - XSO) 

YA  X -GAHX 
H X GAMX 

N X 1 

CALL  FIELOSIHtVA) 

00  10  Jsl,6 
OTNOU)  * (O.Ot  0.0) 

10  QFCURIJ)  X QINTU) 

CALL  FTELDSIYA) 

00  20  Jxl,6 
20  QENOSU)  X QINTIJ) 

CALL  FIELDS  lYS) 

00  30  Jxl,6 

QENOSU)  X QENOSU)  * QINTU) 

30  QTOTALU)  x IQENOSU)  * 4.0»QFCURU)  ) « H / 3.0 
40  00  50  Jxl,6 
50  QYTNTEGU)  < QTOTALU) 
y = H X H / 2.0 
N * 2 • N 

00  60  Jxl.6 

OTNOU)  X QTNOU)  * OFOURU) 

60  QFOURU)  X 10.0,  0.0) 

1 = 0 

70  I = I ♦ 1 

CALL  FIELDS  IY4YA) 

00  80  Jxl,6 

80  OFOURU)  X OFOURU)  ♦ QINTU) 

V * Y ♦ H ♦ H 
IFII  .LT.  N)  GO  TO  70 
IFLAG  X 0 
DO  90  Jxl,6 

QTOTALU)  X IQENOSU)  *■  2.0»QTNOU)  ♦ 4.0*QFOURU)  ) « H / 3.0 
QO<^NOM  X QTOTALU) 

IFICABSIQOENOM)  .LT.  TOL)  QOENOM  x CHPLXITOL.  0.0) 

IFICABSI  lOYINTEGU)  - QTOTALU)  )/OOENOH)  .GT.  TOL)  IFLAG  x i 
90  CONTINUE 

IFIIFLAG  .50.  1)  GO  TO  40 
00  100  Jxl,6 
100  QYTNTEGU)  « QTOTALU) 

RETURN 

END 


Ik 
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SUBROUTINE  FIELDS  IV) 

IMPLICIT  COMPLEXIQ) 

COMHON/ONE/XOt  0NE02F.  ONEO<iF,  ZOt  K*  VO*  OOV2SQ>  QCONST 
C0MH0N/QQ0/QINTI6) 

COHPON/XXX/XOIFF,  XOIFFZ,  XSQ*  X 

REAL  K 

YSO  « V**2 

VOIFF  * VO  - V 

V0IFF2  » V0IFF**2 

TEPMl  a 20  - ONEO«*F*IXSO  ♦ VSO) 

RSO  » XOIFF2  ♦ VOIFF2  ♦ TERM1**2 
R a SORT (RSO) 

ONECR2  a 1.0  / RSQ 
ONFOR  a 1,0  / R 

QCFXP  a CEXPICNPLXIO.Ot  'K^IR  ♦ TERHl  • 20)  ) I 
QPHI  a QCEXP  • CMPLXIONEORt  K)  * 0NC0R2 

OTHETA  a QCEXP  • CMPLX(3.0*ONECR2-K**2»  3. 0*X*ONEOR) •CK£0R2*0NE0R 
TEPM2  a 70  - ONEO<»F  * IXSQ  - V*IV  - 2.0«V0)  ) 

QINTIl)  a QPHI  • TERM2 
QINTI3)  a -QPHI  • XOIFF 
OINTI2)  a -OINTI3)  • V • CNE02F 

QINTI4)  a XOIFF  • IVOIFF  ♦ ONE02Fav*TERNl)  • OTHETA  • QCONST 
OINTI5)  a I2.0*QPHI  - OTHETA* I TERH2*TERN1  ♦ X0IFF2)  )*OCONST 
QINTIE)  a 12. 0*QPHI*V*CNEC2F  - QTHETA* IX0IFF2*V*0NE02F  - VOIFF* 

1 TFRH2)  )*OCONST 
RETURX 
ENO 


